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I. MODEL 

The richness of an electromagnetically active medium is often measured by the variety of 
structures that it can support. Such structures, in turn, are created because of the nonhnear 
response of the medium, for instance, to the impact of a large amphtude electromagnetic 
wave. Naturally the properties of the structure (its shape, its content, its stability, its 
angular momentum etc.) are dictated by the type of nonlinearity that can arise in the 
medium. The discovery or identification of a new nonlinearity type, then, opens up a new 
era of investigation, even, discovery. 

In this paper we work out some of the consequences of a new focusing-defocusing nonlin- 
earity belonging to the general class of saturating nonlinearities (whose magnitude tends 
to a constant as the wave amplitude becomes large). Saturating nonlinearities seem to ap- 
pear, inter alia, in theories of large amplitude wave propagation in pair plasmas (plasmas 
whose main constituents have equal mass and opposite charge in which the pair 

symmetry is broken by some physical mechanism. For instance, a small amount of Bary- 
onic matter (protons) may break the symmetry of an electron-positron (e-p) plasma in the 
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MEV era of the early universe [51, 16|, 17|, 18| . In recently created pair ion (PI) plasmas in the 
laboratory, a variety of symmetry breaking mechanisms like the small contamination by a 
much heavier immobile ion, or a small mass difference between the two constituent species, 
could exist ay , lul, ll2|, ll3|. Asymmetries originating in small temperature differences in 
the constituent species may be always available for structure formation: in the laboratory 
such a temperature difference could be readily engineered (in a controlled way) and there 
are reasons to believe that species temperature difference could exist in cosmic and astro- 
physical setting where one encounters e-p plasmas. It is in this latter setting that a new 



type of nonlinearity 
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appeared while deriving the wave equation (in parabolic approximation) [l| 



describing the nonlinear evolution of the vector potential of an electromagnetic pulse propa- 
gating in an arbitrary pair plasma with temperature asymmetry. Following assumptions and 
notations are necessary in order to put equations ([I]) and ([2]) in perspective: A is the slowly 
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varying amplitude of the circularly polarized EM pulse ~ A (x + y) exp(ikoz — c^ot) with 
mean frequency t^o and mean wave number /cq; V^^ = d"^ /dx^ + d'^ /dy^ is the diffraction 
operator and ^ = z — Vgt is the "comoving" (with group velocity Vg) coordinate. 

Equation ([2]) is written in terms of the dimensionless quantities A = \e\A/{mG{TQ)c'^), 
r = {uJe/c)r, t = uj(,t] where Ue = {'ine^no/my^'^ is the electron Langmuir frequency and 
m is the electron mass. The charges and masses of positive and negative ions are 
assumed to be same (in this paper we mainly concentrate on the specific case of pair plasma 
consisting of electrons and positrons, i.e. q'^ = = = —e~ = \e\ and m+ = m~ = m). 
The equilibrium state of the system is characterized by an overall charge neutrality = 
Uq = no where and Uq are the unperturbed number densities of the positive and negative 
ions respectively. The background temperatures of plasma species are (Tg*" 7^ Tq") 
and mG{z^) = m K-^{z^) / K2{z^) is the "effective mass", [z^ = mc^/T^], where K^, are 
the modified Bessel functions. For the nonrelativistic temperatures (T^ <^ mc^) = 
l + 5T^/2mc^ and for the ultra-relativistic temperatures (T^ ^ rriaC^) G^ = AT^/mc^ ^ 1. 
The smallness parameter e = [G{T^) — G{Tq)]/G(T(I') measures the temperature asymmetry 
of plasma species. For the nonrelativistic temperatures e = 5(Tq^ — TQ)/2mc^ while in 
ultrarelativistic case e = (Tq'" — Tq)/T(^. The numerical factor k = 1/2 for non-relativistic 
temperatures (= 2/3 for ultrarelativistic temperatures). In deriving Eq.([2]) with ([1]), we 
have assumed that the plasma is transparent {ujq ^ 1, Vg ^ 1), and that the longitudinal 
extent of the pulse is much shorter than its transverse dimensions. However, despite of 
dA/d^ ^ V_lA , the second and the third terms in Eq.(l2]) can be comparable due to the 
transparency of the plasma (tuo ^ !)• 

With self-evident renormalization the equation ([2]) can be written as: 

^^+|^+viA+/(iAn.A=o, (3) 

where the nonlinearity function is now following 

which has an unusual feature - in the ultrarelativistic limit (|y4p ^ 1) it tends to be 0. 

Note that the nonlinear refraction index for the considered system can be written as 
8nni = f{I)i where / = |y4p is the intensity of the EM field. The medium is a self-focusing 
{d{Snni)/dI > 0) provided / < 1 while for higher intensities (/ > 1), it becomes defocusing 
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{d(Snni)/dI < 0). For the localized intense EM pulse with the peak intensity > 1 the 
medium becomes defocusing at the peak while remaining focusing at the wings of the EM 
pulse intensity profile. 

In we have demonstrated that Eq.(l3]) supports existence of the stable solitonic struc- 
tures for any spatial dimensions {D = 1,2,3). Such "light-bullets" exist provided that the 
amplitude of the solitons is lower than certain critical values (for instance, in 1-dimensional 
(ID) media < A^cr — 1-4). It is important to emphasize that at —>■ A^cr the profile 
of the central part of the soliton flattens and widens at the top. The existence of fiat-top 
soliton can be explained by the peculiarities of our focusing-defocusing nonlinearity implying 
that the pulse top part with A > 1 entered the defocusing region having the tendency of 
diffraction while the wings of the soliton are in focusing region preventing the total spread 
of the pulse. 



II. FORMATION OF VORTICES 



In this section we examine the possibility of the formation of two-dimensional stable 
soliton-structures carrying a screw type of dislocation, i.e., optical vortices. The generation, 
propagation, and interaction of optical vortices in nonlinear media has been a subject of 
extensive studies (see for review iJ]). In a self-defocusing medium an optical vortex soli- 
ton (OVS) is (2 + l)-dimensional (two transverse dimension and a time) stationary beam 
structure with phase singularity. An OVS is a dark spot, i.e., a zero intensity center sur- 
rounded by a bright infinite background. Self-focusing media also support localized optical 
vortex soliton solutions (LOVS) with phase dislocation surrounded by the bright ring. In 
self-focusing medium LOVS are unstable against symmetry breaking perturbations that lead 
to the breakup of rings into filaments |:15!]. 

Our nonlinearity (jl]) has focusing-defocusing features, hence, one could expect that for- 
mation of both OVS and LOVS solutions is possible in the medium. Such statement can be 
augmented by the results of [l^ where focusing-defocusing model of the media was postu- 
lated to be cubic-quintic medium with sign-changing nonlinearity (/(|Ap) = — \A\^). 
In contrast to the cubic-quintic models the saturation nonlinearity derived by us in jl| is 
not sign-changing - it is the focusing-defocusing one. To verify this expectation we as- 
sume that the pulse is sufficiently long and effects related to the group velocity dispersion 
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(~ cug ^ d'^A/d^^) can be ignored. 

Introducing polar coordinates (r, 9) in {x^y) plane, we look for solutions of Eq.([3]) in the 
form of 

A = A(r) exp(iAt + im6'), (5) 

where integer m defines the topological charge of vortices and A is the nonlinear frequency 
shift. The ansatz ([5]) converts Eq. ([3]) to the ordinary differential equation 

dr'^ r dr (1 + A'^y ~ ' ^ ^ 

We have used numerical methods to find the localized solutions of ([6]). It is possible to 
map the equation in the {A, A^) plane (phase plane) and show that it admits both OVS and 
LOVS solutions. LOVS can exist in the form of infinite number of discrete bound states 
with An{r) (n = 1, 2, ...) where n denotes the finite r zeros of the eigenfunction. 

In what follows we consider only the lowest order (lowest radial eigenmode) solution of 
Eq.([SD {n = 1). For nonzero m (the case we are interested in), the ground state LOVS 
is positive, has a node at the origin r = 0, reaching a maximum, and then monotonically 
decreases with increasing r. Such localized solution exists if A > with the following 
asymptotic behavior: Ar-,o — > r'^'Ao and Ar^^o exp{—r^/X)/y/r , where Aq is a 
constant which measures the slope of A at origin. OVS solutions have the same asymptote 
for r — > 0, while for r — > oo the amplitude has a nonzero value A{r) = Aoo+rn^/ [r"^ f'{Aoo)) 
. Here, A = /(^oo) and OVS exists when /'{Aoo) < , i.e. when background intensity of 
the soliton (far beyond the vortex core) is relativistic A^o > 1. In dimensional units this 
condition corresponds to the negative slope of the nonlinear refractive index {d6nni/dl < 
0), i.e., in the asymptotic region of the solution the medium is defocusing. It is easy to 
demonstrate \v\ that the constant background field with A^o > 1 is modulationally stable. 

A shooting code was used to numerically solve Eq.(l6]). To better understand the results 
of simulations we use the analogy with a nonconservative motion of a particle. For this 
purpose one can rewrite the Eq.([6]) as 

(7) 



d_ 

dr 



''4 I + ViA) 



dr 



r"^ dr r \ dr 



where the "effective potential" is V{A) = -XA^ + ln(l + A'^) - A^/{1 + A^) . The profile 
of the potential V{A) for different values of A is presented in Fig.l. The potential has 
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the maxima at the points A = and A^ax = y 1 — 2A + Vl — 4A /2A . The bounded 
solution is possible only in the case < A < 0.25 while Amax > 1- 

The OVS solutions correspond to a particle beginning its motion at origin {A = 0) with 
certain initial Aq (which can be termed as a velocity (if m = 1) or acceleration (if m = 2 
and so on) and dissipating its initial energy approaching asymptotically potential maximum 
at Amax- Thus, the background intensity of OVS is A^o = Amax, it is always larger than 
unity and can become arbitrarily large for A — 0. We also found out that OVS solutions 
exist even for 0.25 > A > Acr — 0.2162, i.e., when V{A,^s,x) < (see curve "a" in Fig.l). 
In other words, the effective particle cannot overpass but only approach asymptotically the 
lower potential maximum. 

The numerical solutions of the nonlinear equation ([6]) for m = 1, 2 and 3 are shown in 
Fig. 2. As expected, the soliton-like solutions evidently go to zero as r"^ for small r, and 
reach an m- independent asymptotic value predicted above. In Fig. 3 curve "a" displays the 
dependence of the field derivative at the origin, Ao , as a function of the nonlinear frequency 
shift A for m = 1 case. One can see that Aq is a growing function for small A-s . 
For small A-s the position of the potential maximum "moves" to larger values of A and, 
consequently, "particle" needs to have larger initial "velocity" [Aq) to reach the maximum. 

In contrast to OVS, the LOVS solutions correspond to the particle returning back asymp- 
totically to the initial position at A = 0. It seems obvious that due to the "frictional" 
motion particle can not make its way back if A > A^- Thus, LOVS may exist in the range 
< A < Acr while its amplitude (in contrast to OVS) is a growing function of A. Such 
dependence is obtained numerically and is presented in Fig. 4. for single-charged vortices 
(m = 1). One can see that the amplitude of the LOVS (A^) is bounded from above by 
certain critical value for Acr{— 1-5). Thus, in contrast to OVS the localized vortex can be 
just moderately relativistic. Notice that for 0.16 < A < Acr the amplitude of the LOVS 
(Am) varies in the range 1 < Am < Acr ■ For the top part of such solution (with A{r) > 1) 
the medium is defocusing while remains focusing for lower intensity wings of the structure. 
Consequences of this fact can be seen in Fig. 5 where profile of LOVS is given for variety 
of A. With increase of A the central part of the LOVS flattens and widens converging to 
the OVS. In principle, it is possible to create flat-top LOVS with a large transverse width. 
Convergence of LOVS to OVS can be better seen in Fig. 3 where the curve "b" corresponding 
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to LOVS almost coincides with curve "a" near the point A ~ Acr- Similar behavior of the 
solutions can be obtained for the vortices with higher charge (m = 2,3,..), however, the 
corresponding figures we do not present here for brevity. 

III. STABILITY OF SOLUTIONS 

Are these solitonlike solutions stable? 

The intensity dependent switching from the focusing to defocusing regime can have an 
interesting consequence for the stability properties of the solutions. As it is well established 
[l^ OVSs with m = 1 are stable whereas vortices with a larger value of m may decay into 
the single-charged ones in self-defocusing media. In our case the bulk of the OVS is always 
in the defocusing regime and as we mentioned above the background field is always stable. 
However, near to vortex core the medium becomes focusing. Thus, stability of the OVS in 
our case can not be granted. 

We performed stability analysis solving numerically Eq.(l3]); while simulation (for various 
A-s), initial stationary OVS state was perturbed radially and azimuthally by Gaussian noise. 
Typical picture of the evolution is plotted is Fig. 6. We see that perturbations are quickly 
radiated away and the initial state relaxes to the ground state solution implying that the 
OVS is stable in the whole examined range. 

To verify stability of LOVS we first performed a linear stability analysis. To do so we 
follow procedure developed by [17|] and consider perturbation acting along a ring of mean 
radius r^, , where v4(r^,) = ■ Assuming constant intensity and spatial uniformity for this 
ring, one can rewrite the diffraction operator in ([3]) as V^^ = r~'^d'^ jdfP' and introducing 
azimuthal perturbation with a phase factor \E' = f2t + MQ (where M is an integer) for the 
growth rate of instability we get: 

M [2(1 -AM M^' 

Im(^]) = — Re ; - 

One can see from ([S]) that large amplitude LOVS with Am > 1 is always stable. For the 
lower amplitude case LOVS should decay into Mmax multiple filaments, where Mmax is an 
integer close to the number for which maximal growth rate is maximal. 

In Fig. 7 we plot Im(fi) versus M for A = 0.1 and for m = 1,2,3 . Corresponding 
Am are respectively 0.66; 0.65; 0.63 and r^=6.3; 11.6; 16.9. One should expect that 
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instability will result in breaking of the pulse into the filaments (fragments) with number 
of filaments being respectively 2, 4, and 5 (or 6) for m = 1,2,3. These filaments have to 
conserve the total angular momentum. Since the fusion of filaments is not possible due to 
the topological reasons, they can eventually spiral about each other or fly off tangentially to 
;he initial ring generating bright solitonic structures found for index saturation nonlinearity 



Our numerical simulations for A.^ < 1 give evidence of a quickly developing instability 
in agreement with predictions of linear stability analysis. Indeed, in Fig. 8 we present the 
results of simulations for the LOVS with m = 1 and 2. The LOVS with m = 1 breaks into 
two filaments while for m = 2 the breaking into 4 filaments takes place. The filaments 
are running away tangentially without spiraling. All filaments like spatial solitons remain 
stable. Most interesting is the situation when amplitude of LOVS is larger than unity. 

In Fig. 9 we present the evolution of LOVS both for m = 1 and m = 2 ; corresponding 
amplitudes for the soliton solutions are (A = 0.2) : = 1-39 and A^ = 1.37 , 

respectively. The initial input LOVS solution was modulated by a Gaussian noise (the level 
of noise was 5%). One can see that breaking of the LOVS does not take place. In order to 
be sure that some very slow instability is not developing the simulations were carried out 
until t = 4000, i.e., for 130 soliton period Tgoi = 27r/A ~ 30. Thus, single and multi-charged 
LOVS become stable for large amplitudes. 

At the end we would like to emphasize that if for single-charged LOVS we are confi- 
dent about its stability for multi-charge LOVS (as well as for OVS) one should be careful. 
Indeed, from general topological reasons the multi-charged vortices are supposed to be un- 
stable and they should break into single charge vortices. However, we could not observe 
in our simulations such breaking. It is possible that this instability is very slow (sub- 
exponential/algebraic) and, as a result we obtained that multi-charged vortices are very 
long-lived objects and practically stable. 

In our consideration the effects related to the group velocity dispersion and corresponding 
reshaping of the radiation have been ignored. However, one can generalize our results by 
keeping the term ~ d'^A/d^'^ in Eq.(l3]). In transparent plasma case this term can affect the 
long time dynamics of self-guiding vortex solitons. In particular due to weak modulation 
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instability [18] the self-trapped beam eventually will break into a train of spatiotemporal 
solitons, i.e. the "light bullets". However, due to the topological reasons the vortex lines 
should survive the structural changes. We expect that instability will result in generation of 
fully localized bullets of vortex solitons (the spinning-bullets). Dynamics of formation and 
stability of such structures is beyond of the scope of current paper. 

IV. CONCLUSIONS 

The asymmetries originating in small temperature differences in the constituent species 
of electromagnetically active medium may be always available for structure formation both 
in laboratory and cosmic/astrophysical settings. In present paper we have shown that this 
asymmetry defines the specific properties of the structure due to the new type of the non- 



linearity that can arise in the medium (derived in our earlier paper [l]). We found that the 
pair plasmas of any dimension with "asymmetry" in initial temperatures can support the 
stable large amplitude optical vortex and localized vortex solitons. Localized structures for 
certain parameters may have the flat-top shapes. The coexistence of LOVS and OVS solu- 
tions and their stability in such medium is due to the specific form of saturating nonlinearity 
switching from the self-focusing to the self-defocusing regime and vice versa. This fact is 
rather interesting specifically for laboratory conditions. 
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FIG. 1: The "effective potential" versus the amphtude for different values of the nonlinear 
frequency shift A. The curve "a" corresponds to A > \cr — 0.2162, the curve "b" has A = Xcr, and 
for the curve "c" < A < Xnr ■ 




FIG. 2: Profiles of OVS-s for m = 1, m = 2, m = 3; nonlinear frequency shift A = 0.16. 
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FIG. 3: Aq versus A for m = 1; curve "a" corresponds to OVS while curve "b" - to LOVS. 




FIG. 4: The effective eigenvalue A versus soliton amplitude for m = 1. 




FIG. 5: Profiles of soliton solutions. Curves "a", "b" and "c" correspond to LOVS with A 
0.005; 0.16; 0.205, respectively. Curve "d" corresponds to OVS for A = 0.205. 
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FIG. 6: The dynamics of initially perturbed OVS-s: plots are chosen for different time-moments 
t = 0; 100; 1000 . 
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FIG. 7: Instability growth rate Im{Q) versus M for A = 0.1 for different topological charges m. 
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FIG. 8: Vortex dynamics (for different time-moments) when A = 0.1: the left panel - for m = 
1) Ajnax = 0.66 , the vortex splits into 2 filaments; the right panel - for m = 2, Ajnax = 0.6580, 
the vortex splits into 4 filaments; the filaments are running away tangentially. 
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FIG. 9: Vortex dynamics (for different time-moments) when A = 0.2, the vortex is robust towards 
perturbations; the left panel is for m = 1, A^ax = 1-386 while the right panel is for m = 2, = 
1.3729. 



